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Abstract 

We consider globally regular and black holes solutions for the Einstein- Yang-Mills system with neg- 
ative cosmological constant in d— spacetime dimensions. We find that the ADM mass of the spherically 
symmetric solutions generically diverges for d > 4. Solutions with finite mass are found by considering 
corrections to the YM Lagrangean consisting in higher therm of the Yang-Mills hierarchy. Such systems 
can occur in the low energy effective action of string theory. A discussion of the main properties of the 
solutions and the differences with respect to the four dimensional case is presented. The mass of these 
configurations is computed by using a counterterm method. 

1 Introduction 

If we allow for a negative cosmological constant, the solution of the matter free Einstein equations possessing 
the maximal number of symmetries is the anti-de Sitter (AdS) spacetime. Recently a considerable amount 
of interest has been focused on solutions of Einstein equations with this type of asymptotics. This interest 
is mainly motivated by the proposed correspondence between physical effects associated with gravitating 
fields propagating in AdS spacetime and those of a conformal field theory on the boundary of AdS spacetime 

mm 

In view of these developments, an examination of the classical solutions of gravitating fields in asymp- 
totically AdS (AAdS) spacetimes seems appropriate. Recently, some authors have discussed the properties 
of gravitating SU(2) nonabelian fields with a negative cosmological constant A Considering the case 

of four spacetime dimensions, they obtained some surprising results, which are strikingly different from the 
results found in the asymptotically flat case. For example there are solutions for continuous intervals of the 
parameter space, rather then discrete points. The asymptotic values of the gauge potentials are arbitrary 
and there exist solutions supporting magnetic and electric fluxes without a Higgs field. Some of these solu- 
tions are stable against spherically symmetric linear perturbations. The literature on AAdS solutions with 
nonabelian fields is growing steadily, including stability analyses jS] |S| , the study of configurations with a 
NUT charge [7] , topological black holes with nonabelian hair [5] as well as axially symmetric generalisations 
[3 E|- The existence of these solutions invalidates the AdS4 version of the no hair conjecture, which states 
that the black holes are completely characterised by their mass, charge and angular momentum. 

However, all these studies approach the case of a four dimensional AAdS spacetime, and relatively little 
is known about higher dimensional AAdS solutions with non Abelian matter fields. Practically all that is 
known for d > 4 is the five dimensional nonabelian SU(2) solutions discussed in At the same time 

gauged supergravity theories playing an important role in AdS/CFT, generically contain non Abelian matter 
fields in the bulk, although in the literature only Abelian truncations are considered, to date. Thus, the 
examination of higher dimensional gravitating non Abelian solutions with A < is a pertinent task. 

Higher dimensional asymptotically flat solutions of the Einstein- Yang-Mills (EYM) equations have re- 
cently been the subject of several studies. As found in |12) . in five spacetime dimensions, the particle 
spectrum obtained by uplifting the d = 4 flat space YM instantons become completely destroyed by gravity, 
as a result of their scaling behaviour. However, by adding higher order 1 terms in the Yang-Mills (YM) 
hierarchy this obstacle due to scaling is removed and solutions in higher dimensions can be found. Such 
regular, static and spherically symmetric solutions in spacetime dimensions d < 8 were presented in |13| . 

1 The only higher order curvature terms considered in this paper are those constructed from a 2p-form field such that the 
Lagrangean contains velocity square fields only. 



and for d = 5, both globally regular and black hole solutions of the EYM system were found in |14| . The 
properties of these solutions are rather different from the familiar Bartnik-McKinnon (BK) solutions ^j] to 
EYM in d = 4, and are somewhat more akin to the gravitating monopole solutions to EYM-Higgs (EYMH) 
|17| . This is because like in the latter case |TJ|, where the VEV of the Higgs field features as an additional 
dimensional constant, here J3| also additional dimensional constants enter with each higher order YM cur- 
vature term 2 . They are however quite distinctive in their critical behaviour. The typical critical features 
discovered in [131 114j have recently been analysed and explained in |19| . 

These results can be systematically extended to all d > 5 and one finds that no finite mass solutions can 
exist in EYM theory unless one modifies the non Abelian action by adding higher order curvature terms 3 
in the YM 2p-form curvature F(2p), the p = 1 case being the usual 2-form YM curvature. Without these 
higher order YM terms, only vortex-type finite energy solutions |12) exist, describing effective systems in 3 
spacelike dimensions, and with a number of codimensions. 

It is the purpose of this paper to examine the corresponding situation in higher dimensions, in the 
presence of a negative cosmological constant. It is interesting to inquire whether the introduction of a 
negative cosmological constant to these higher dimensional EYM models will lead to some new effects as it 
does in the d = 4 case, due to the different asymptotic structure of the spacetime. In the first place this 
would lead to our understanding of how the behaviour of EYM theory depends on the dimensionality of the 
spacetime. But such higher dimensional AAdS solutions might be relevant to superstring theory, namely 
in the context of solutions to various supergravities containing nonabelian matter fields. Here however we 
restrict our considerations to the simplest case of systems consisting only of gravitational and YM fields, 
namely to higher dimensional EYM model. In particular, we restrict to Einstein- Hilbcrt gravity and the 
first two members of the YM hierarchy, and hence to d < 8, augmented with negative cosmological constant 
in d spacetime dimensions. As in [312], f° r the d = 4 case, we seek static spherically symmetric solutions 
in the d — 1 spacelike dimensions. We find both globally regular and black hole solutions with finite ADM 
mass. Unlike in the d = 4 case however, we find that for d > 4 and a negative cosmological constant, the 
properties of the AAdS solutions do not differ qualitatively from the asymptotically flat case. 

Our strategy is to first consider the usual YM model, namely the p = 1 member of the YM hierarchy, 
or the square of the 2-form curvature F(2). We present an argument for the absence of solutions with 
reasonable asymptotics for any spacetime dimension d > 5. Although the EYM equations in this case 
present solutions approaching asymptotically the AdS background, the mass generically diverges. This can 
be seen as a simple version of the no hair theorem, holding for the EYM system in d > 4 dimensions. In 
other words, Schwarzschild-AdS black hole is the only static, spherically symmetric solution of the EYM 
system with finite mass. This is presented in Section 2, where in addition we have considered the special 
case d = 3, extending Deser's analysis [2]] for A < 0. 

In Section 3 we introduce the higher order YM hierarchy models featuring the terms F(2p), p > 2. We 
derive the classical equations subject to our spherically symmetric ansatz, and present a detailed numerical 
study of of both regular and black hole solutions. As is the case with the usual EYM system, where the 
existence of finite mass regular solutions leads to the violation of the no hair conjecture, here too, these 
solutions are explicit counterexamples to this conjecture in AdS^ spacetime. 

One may ask about the possible relevance of these higher dimensional configurations within the AdS/CFT 
correspondence. In Section 4 we compute the boundary stress tensor and the mass and action of the solutions 
in a number of spacetime dimensions up to eight. In five and seven dimensions, the counterterm prescription 
of [21] gives an additional vacuum (Casimir) energy, which agrees with that found in the context of AdS/CFT 
correspondence. A counterterm based proposal to remove the divergences of a F(2) theory, such that the 
mass and action be finite, is also presented in Appendix B. We give our conclusions and remarks in the final 
section. 

2 Like the gravitating monopoles, these have gravity decoupling limits except in d = 5 (and in d = 4p + 1 modulo 4), and in 
all odd spacetime dimensions, the flat space solutions are stabilised by a Pontryagin charge analogous to the magnetic charge 
of the monopole, provided that the representations of the gauge group are chosen suitably. In all even d however, they are like 
the BK solution in that they are not stabilised by a topological charge and are likewise sphalerons 115) 

' j in principle higher order terms with the desired scaling can be chosen to consist both of the YM and the Riemann curvatures, 
but in practice we restrict to the YM hierarchy. The reason will be explained in section 2.5. Besides, it was found in 1131 that 
the inclusion of Gauss-Bonnet terms does not result in any new qualitative features to the solutions. 
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Everewhere in this paper we employ the notations and conventions of |13) . 



2 The F{2) model 

2.1 The action principle 

We start with the following action principle in d— spacetime dimensions 

M d ^ (lab** - 2A ) + c -) ~ L d3xV ^ hK > « 

where R is the Ricci scalar associated with the spacetime metric A = — (d — l)(d — 2)/(2l 2 ) is the 
cosmological constant and G is the gravitational constant (following we define also n — 1/(8ttG)). 
The matter term in the above relation 

C m = ~n tr F^F^ (2) 

is the usual F 2 nonabelian action density, Fp V — d^A v — d„A^ — i[Afj,,A v ] being the gauge field strength 
tensor. Here t± is the coupling constant of the model (in the usual notation n = \/g 2 ). 

The last term in is the Hawking-Gibbons surface term JSji where K is the trace of the extrinsic 
curvature for the boundary dM. and h is the induced metric of the boundary. Of course this term does not 
affect the equations of motion but it is relevant for the discussion of the mass and the action of the solutions, 
in Section 4. 

The field equations are obtained by varying the action Q with respect to the field variables , A^ 

Rfiu - 7j9 t iiyR + j^-g^iy = 87rGT Ml/ , (3) 
V^-iiA^F^] = 0, 

where the energy momentum tensor is defined by 

T^ v = tr F^F^g " 3 - ^tr F af3 F af3 . (4) 

2.2 The general Ansatz 

For the case of a G?-dimensional spacetime, we restrict to static fields that are spherically symmetric in the 
d — 1 spacelike dimensions, with a metric ansatz in terms of Schwarzschild coordinates 

ds2 = W) + rHv?d - 2 - ^M^M^ 2 ' ( 5 ) 

with c?ri c j_ 2 the d — 2 dimensional angular volume element and 

2m(r) r 2 , . 

the function m(r) being related to the local mass-energy density up to some d— dependent factor. 

As discussed in ^21> the choice of gauge group compatible with the symmetries of the line element ||5) 
is somewhat flexible. In J3| the gauge group chosen was SO(d), in d dimensions. But the gauge field of 
the static solutions in question took their values in SO(d— 1). Thus in effect, it is possible to choose SO(d) 
in the first place. Now for even d, it is convenient to choose SO(d) since we can then avail of the chiral 
representations of the latter, although this is by no means obligatory. Adopting this criterion, namely to 
employ chiral representations, also for odd d, it is convenient to choose the gauge group to be SO(d — 1). 
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We shall therefore denote our representation matrices by SO± (d) , where d = d and d = d — 1 for even and 
odd d respectively. 

In this unified notation (for odd and even d), the spherically symmetric Ansatz for the SO± (devalued 
gauge fields then reads 

^o = 0, A, = (i^) Elf i, , ^=-l(}^)[T i ,T j] . (7) 

The T's denote the <i-dimensional gamma matrices and 1, j = 1, 2, d — 1 for both cases. 
Inserting this ansatz into the action QJ, the EYM field equations reduce to 

= (r d - 4 <jNw')' -{d-3)r d - & a{w 2 - l)w, (8) 

- - T \^(N W - + { d- 3)^!), (9) 

(10) 



fj Tl W 

which can also be derived from the reduced action 



2 V 2r 2 

/2 



5-/*. (m' - ^ 4 (iW 2 + (d- 3)^3^)) • (11) 

For a -F 2 theory, the constants k, ti can always be absorbed by rescaling r — * cr, A — * A/c 2 and m — > mnc d ^ 3 , 
with c = y ti/k. 

The above differential equations have two analytic solutions. One of them with 

w(r) = ±1, m(r) = M, cr(r) = 1 (12) 

corresponds to Schwarzschild-AdS spacetime. For iu(r) = we find a non Abelian generalisation of the 
magnetic- Reissner-Nordstr0m solution with a(r) — 1 and 

m(r) = M + — logr if d = 5, and m(r) = M + n } d ~ 3 ) r ^~5 for d ^ 5i ( 13 ) 
2 4(a — 5) 

Mo being an arbitrary constant. We can see that, although these solutions are asymptotically AdS, the mass 
defined in the usual way diverges. 



2.3 d = 3 

The (2+1) dimensional case is rather special. Three dimensional gravity has provided many important clues 
about higher dimensional physics. This theory with A < has non-trivial solutions, such as the BTZ black- 
hole spacetime |22| . which provide an important testing ground for quantum gravity and the AdS/CFT 
correspondence. Many other types of 3c? regular and black hole solutions with a negative cosmological 
constant have also been found by coupling matter fields to gravity in different ways. 

However, as proven in [201, there are no d = 3 asymptotically flat static solutions of the EYM equations. 
The arguments in can easily be generalised for the AAdS case. We notice that for d = 3, the YM 
equation (jSJl implies the existence of a first integral w' = ar/(aN) with a an arbitrary real constant. 
Therefore, assuming AdS3 asymptotics, w' decays asymptotically as 1/r which from JJjJ implies a divergent 
value of m(r) as r — > oo. However, similar to the A = case |23| . this argument does not exclude the 
existence of nontrivial solutions of the field equations. 

Here we should remark that since for d = 3we are dealing with SO(d — 1) = SO (2), the gauge group is 
Abelian and we recover Einstein-Maxwell theory with a negative cosmological constant, whose solutions are 
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known in the literature. The corresponding solution with a vanishing electric field was found by Hirschmann 
and Welch [21] and has a line element 4 



r 



.2 



dr 2 



ds 2 = 



2 + c 2 \og\r 2 /l 2 



M\ (r 2 /£ 2 - M) 



+ (r 2 + c 2 log |r 2 /^ 2 - M\)dip 2 - {r 2 /£ 2 - M)dt 2 , 



(14) 



r 



the magnetic potential being 



w(r) =w + —= log |r 2 /^ 2 - M|, 



(15) 



with wq , c and M arbitrary real constants, the BTZ metric being recovered for c = (see also [251 126) for 
more details on this solution). 

One can see that although the quasilocal mass defined in the usual way diverges as r — > oo, the metric 
still approaches the AdS3 background. However, a similar problem appears for other d = 3 AdS solution 
e.g. for the electrically charged BTZ black hole [22], or for a self- interacting scalar field minimally coupled 
to gravity [^ , in which cases it was possible to find a suitable mass definition. We expect the formalism 
developed in those cases to work also for the Hirschmann- Welch solution l|14|) - (|15fl . but this lies outside the 
scope of the present work. 



Four dimensional black hole solutions of the equations (|5|)- (|10|l have been found in [3J, the globally regular 
counterparts being discussed in 0]. Differing from the asymptotically flat case, for A < there is a continuum 
of solutions in terms of the adjustable shooting parameter that specifies the initial conditions at the origin 
or at the event horizon. As a new feature, the asymptotic value of the gauge function wq is arbitrary. The 
spectrum has a finite number of continuous branches, depending on the value of A. When the parameter 
A approaches zero, an already-existing branch of solutions collapses to a single point in the moduli space. 
At the same time new branches of solutions emerge. A fractal structure in the moduli space has been 
noticed [2U- There are also nontrivial solutions stable against spherically symmetric linear perturbations, 
corresponding to stable configurations. The solutions are classified by non-Abelian magnetic charge and the 
ADM mass. 

Note also that the d = 4 EYM solutions with a negative cosmological constant A = — 3/ti have some 
relevance in AdS/CFT context. As proven in [25], for this value of the cosmological constant, an arbitrary 
solution {g^v, A^) of the four dimensional EYM equations gives a solution of the equations of motion of 
the d = 11 supergravity. Based on this observation, an exact BPS-type EYM solution has been constructed 
in |1U| . However, similar to some supersymmetric solutions in Einstein-Maxwell theory with A < 0, this 
A = — 3/ti configuration presents a naked central singularity. 



As discovered by Coleman j3U] and Deser [3T], there are no flat space static solutions of the YM equations, 
except for d = 5. However, the inclusion of gravity may change this picture, as seen from the famous d = 4 
asymptotically flat Bartnik-McKinnon solutions ^BJ. In this case, the repulsive YM force is compensated 
by the attractive character of the gravity, and as a result we find both regular and black hole unstable 
configurations (see [32] for a fairly recent survey). As found in ^2] the d = 5 particle- like solutions are 
destroyed by gravity, their mass diverging logarithmically, while w(r) presents an infinite number of nodes. 
The AAdS couterparts of the d — 5 asymptotically flat solutions are discussed in . Although approaching 
asymptotically the AdS$ background, the mass of these configurations also diverges logarithmically. 

As conjectured by several authors, this result extends to higher dimensions. Following the approach in 
we prove in Appendix A the nonexistence of asymptotically flat or AdS solutions with a finite mass 

4 One can also solve directly the field equations ISI- 11UI . but the solution takes a much more complicated form for the metric 
ansatz 



2.4 d = A 



2.5 d>5 
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in a F(2) EYM model given by @ for any spacetime dimension d > 5 (see also the discussion in Section 
3.3). Therefore, the d— dimensional Schwarzschild-AdS configuration is the only finite mass solution of the 
equations JSJ-^UJ) and a simple version of the no hair theorem seems to hold for the F(2) (usual) EYM 
system in d > 4. 

We should remark that in deriving this result we assumed implicitely that m(r),a(r),w(r) are smooth 
functions approaching finite values as r — > oo. A divergent asymptotic value of m(r) invalidates the proof 
presented in Appendix A and also the virial argumnets in Section 3.3. Therefore we cannot exclude the 
existence of spherically symmetric, nontrivial solutions of the field equations for any d > 4. However, the 
mass of these solutions generically diverges, although the spacetimes are still AAdS. The work of Ref. JJ| 
presents an extensive discussion of such AAdS solutions for d = 5. Both regular and black hole solutions 
exist in d = 5 for compact intervals of the parameter that specifies the initial conditions at the origin or at 
the event horizon. Differing from the A = case, the gauge field function w(r) does not oscillate between 
1 and —1 but approaches asymptotically some finite value wq, the node number being finite. The masses 
of these solutions behave asymptotically as (wq — l)logr, with all wq — ±1 solutions corresponding to pure 
gauge configurations. 

The results we found by solving numerically the equations l(8*)l-(|10[l. for d = 6,7,8 and several negative 
values of A confirm that this is a generic behaviour for d > 4. The corresponding boundary conditions at 
the origin (or event horizon) are found by taking P = 1 in relations l|26(l . Q27JI given in Section 3. Except for 
a divergent value of m{r) as r — ► oo, according to 

m(r) = Mo + - 1) V- 5 , (16) 

the properties of these solutions are very similar to the the more familiar d = 4 case. For d > 5, the 
asymptotic value wq of the gauge field function w is also arbitrary, being fixed by the initial parameters 
u>"(0) or w(rh) respectively, wo = ±1 corresponding to pure gauge configurations. Solutions for a compact 
interval of these parameters were found to exist, the general structure being A-dependent. Solutions with 
nodes in w(r) were also found. Typical d — 6 configurations with a regular origin are presented in Figure 
1, for A = —0.01 and three different values of b = —w"(0)/2. One can see that the mass function diverges 
linearly while cr(r) and w(r) asymptotically approach some finite values. In Figure 2 we plot the parameters 
M (appearing in (|16|) . which in Section 4.3 we argue that it can be taken as the renormalised mass of the 
solutions), wo, the value ao of the metric function a at the origin and the minimal value N m of the metric 
function iV as a function of b for a family of d — 6 AAdS solutions with A = — 1. This branch ends for some 
finite value of b, where <r(0) — > 0. Black hole solutions have been found as well, presenting the same general 
features. Here also we find a continuum of solutions with arbitrary values of Wq, the relevant parameter 
being the value of the gauge potential at the event horizon. Similar to d = 4, solutions appear to exist for 
any value of the event horizon radius. 

The drawback of the solutions in d > 5 described above is that their ADM masses are divergent, making 
their physical significance obscure (see, however, the discussion in Appendix B). 

One may hope to find a different picture by including some other matter field in the action Q • Such fields 
should interact with the YM sector so as to compensate for the scaling behaviour of the non Abelian fields. 
This excludes the dilaton field, as it can be proven that the latter does not change this nonexistence result. 
Note that d = 5 finite mass spherically symmetric gravitating non Abelian solutions with a Liouville-type 
dilaton potential are known to exist |35| : however these solutions are asymptotically neither flat nor AdS. 

In the next Section, we will remedy this problem of nonexistence by adding higher order YM curvature 
terms F(2p), (p > 1), to (J2J. The main role of these terms is to alter the scaling properties of the action 
density in . But logically, such a role can be played also by altering the gravitational part of QJ , through 
the introduction of higher order (in the Riemann curvature) terms R{p), (p > 1). R(p) here denotes the 
generalised Ricci scalar constructed from the 2p-form antisymmetrised p-fold product of Riemann tensor 
2-forms, the p = 1 member of which gives the Einstcin-Hilbcrt action, the p = 2, the Gauss-Bonnet, etc. 
Here we have eschewed the possibility of employing additional R(p) terms instead of, or together with F(2p) 
YM terms, because in the present work we exclude the participation of fields other than gravitational and 
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Yang-Mills. In particular, the exclusion of the dilaton renders the usefulness of higher order (Gauss-Bonnet) 
gravities trivial from a practical point of view. 

To see this, consider the situation where a term scaling as L~ 2p is needed, i.e. that the term 

must be added to the density in {Tjl. Using the metric Ansatz (|5j), and discarding purely boundary terms in 
the residual one dimensional (spherically symmetric) Lagrangean, the term C( p ) to be added is 



Jr d-2v-2 



r4(l - N) p + (d-2p- 1)(1 - N) p 
ar 



(17) 



in spacetime dimension d, which vanishes identically in dimensions up to d = 2p. This is because in the given 
dimensions the p-th member of the gravitational hierarchy becomes a total divergence. But for given d it is 
necessary to have 2p > d, whence ()17|l trivialises. The situation would be quite different if the dilaton field 
were included together with the higher order gravitational terms, preventing the p-gravity density becoming 
a total divergence. But this option is excluded here and we opt to the exclusive use of higher order p-YM 
densities. 

Yet other fields occurring in supergravities might be considered, e.g. Kalb-Ramond, or totally antisym- 
metric tensor fields. But these being Abelian, their effect would be felt only in given dimensions, or, subject 
to much less stringent symmetries than the spherical. This also is not a flexible option, so we restrict our 
attention to p-YM densities only. 



3 Nonabelian hierarchies 

3.1 The Lagrangean and field equations 

Since finite mass spherically symmetric solutions play a central role in AdS/CFT, it is desirable that the 
nonexistence result presented in Section 2.5 be circumvented. 

A simple way to circumvent these arguments and to find nontrivial solutions is to to modify the matter 
Lagrangean by adding higher order terms in the YM hierarchy. As noted in footnote 1 , these are constructed 
exclusively from YM curvature 2p-forms. For A = 0, asymptotically flat, finite energy solutions of this 
modified EYM system are constructed in [131 114) . 

Such terms as we propose to add are predicted by string theory, and hence provide a link with the 
AdS/CFT correspondence too. But here we are guided predominantly by symmetry considerations and do 
not claim to be employing terms strictly following from superstring theory. The situation concerning higher 
order YM curvature terms in the string theory effective action is complex and as yet not fully resolved. 
While YM terms up to F A arise from (the non Abelian version of) the Born-Infeld action [HE], it appears 
that this approach does not yield all the F 6 terms [HZ] ■ Terms of order F 6 and higher can also be obtained 
by employing the constraints of (maximal) supersymmetry [38j . The results of the various approaches are 
not identical. 

The definition we use for superposed YM hierarchy is 

p i 

£ ™ = - E 2(2rf! Tp ^ Tr F{2P)2 ' (18) 
where F(2p) is the 2p-form p-fold totally antisymmetrised product of the SO(d) YM curvature 2-form F(2) 

Even though the 2p-form l|19|l is dual to a total divergence, namely the divergence of the corresponding 
Chern-Simons form, the density l|18fl is never a total divergence since it is the square of one. But the 2p-form 
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(|19fl vanishes by (anti)symmetry for d < 2p so that the upper limit in the summation in (|18|l is P = i for 
even d and P = ==i for odd d. 

We define the p-stress tensor pertaining to each term in (|18fl as 

T« = Tr F(2p) AlAlA2 ... A2p _ 1 ^(2p)/ lA - A ^- 1 - ^ 5ai „ Tr F(2p) Al Aa ... A2p F(2p) Al A - A2 » . (20) 

For the particular spherically symmetric ansatz considered in Section 2, we express the reduced YM 
Lagrangean arising from i|18[) as 

p 

C - - V J (p) 

>~m — 2_^i YM 

p=l 



(21) 



(p) 

with iy^f given by 



T (P) - r d-2 T r> 
YM ~ r 



(d-2)! 



2-(2p)! (d-[2p+l})\ 
having used the shorthand notation 

W = 

For this general ansatz, we find the field equations 



wp- 1 



(2p)N 



1 dw 
r dr 



+ (d-[2p+l})W 



2 2^ of 



„d-2 



(d-3)! 



£j 2(2p)! (d-[2p+l]) 



■VFP- 1 



(2p)JV 



w 2 - 1 



1 diu 
r dr 



(22) 



(23) 



(d- [2p+l])W 



a 1 (d-3)! 



P 



P-2 



^(d-[2p + l])!(2p-l)f 



(24) 



^dr\(d-[2 P +l})\(2p-l)\ ) ^ (d-[2p+l])!(2p-l 



= 1 " - - p= 

which can also be derived from the reduced action 



. w' 2 d - 2p - 1 



S= dr 



,d-2 



(d-3) 



£j 2(2p)! (d-[2p+l])l 



(2p)7V ( 



ldtoV 
r dr J 



+ (d-[2p+l])W 



(25) 



For T, = 0, i > 2 the equations l|8")l- H10(l are recovered. Note also the a equation decouples and can be treated 
separately. 



3.2 Boundary conditions 

The asymptotic solutions to these equations can be systematically constructed in both regions, near the 
origin (or event horizon) and for r>l. 
The corresponding expansion as r — > is 



w(r) = l-6r 2 +0(r 4 ), 



m(r) = 
a(r) = 



,(d-3)!(46 2 )f 



§2(2p)!(d-[2p+l])! 



„d-l 



0(r 



d+l\ 



(26) 



ctq 



(d-3)! 
2k 



t p (46 : 



(§(2p-l)!(d-[2p+l])! 



)r 2 +0(r 4 ), 
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and contains one essential parameter b (the value of er(0) = a can be fixed by rescaling the time coordinate). 

For black hole configurations with a regular, nonextremal event horizon at r = r^, the expression near 
the event horizon is 

m(r) = m h + m'(r h )(r - r h ) + 0(r - r h ) 2 , 

w(r) = w h +w'{r h )(r-r h ) + 0{r-r h ) 2 , (27) 
cr(r) = a h + <r' h (r - r h ) + 0(r - r h ) 2 , 



where 



m{r h ) = -nr d h ~\l + p, W h = 

m / fr ,l-V ^- 2 (ri-3)! w p N , d-3 (d-l)r h 2m>(r h ) 
m ^-Z. 2 (2p)!(d-2p-2)!^' ^" r„ + ^ 2 K ^" 3 ' (28) 



, = a h (d-3)\ , 2 >A TpWl 1 , = f W h (w 2 - 1) £p=l (2p-l)!(d-2p-2)! 

a ' 1 /sr/, h (2p- l)!(d- [2p + l])!' h N' h r\ v p ^w^ 1 



= 1 (2p-l)!(d-2p-l)! 

the value of the gauge field on the event horizon being the esential parameter. Here the obvious condition 
N'(rh) > imposes some limits on the event horizon radius as a function of Wh for given (ji, A). 

Since the field equations are invariant under w — » — w, one can take w(0) = 1 and w(rh) > without 
any loss of generality. 

For r> 1 we find for both regular and black hole solutions 

m W = ~ 8£ 2 + (29) 

^(d-3) 2 n l 

CTW = 1_ 2«(d-2) ^ + --" 

These boundary conditions are also shared by the asymptotically flat solutions (with a different decay of the 
mass function m(r), however), w = ±1 being again the only allowed values of the gauge function as r — > oo. 
Therefore, we expect to find a qualitatively similar picture in both cases. We will find in Section 5 that the 
constant M in the above relations is the ADM mass up to a d-dependent factor. However, in the discussion 
of numerical solutions we will refer to M to as the mass of the solutions. 



3.3 Further relations 

The form (|25|) of the reduced action allow to derive an useful virial relation. To this end, we use the scaling 
technique proposed in |331 I34| for the case of spherically symmetric gravitating systems. Let us assume the 
existence of a globally regular solution m(r),a(r),w(r) of the field equations with suitable boundary 
conditions at the origin and at infinity. Then each member of the 1-paramctcr family 

m\(r) = m(Ar), a\(r) = a(\r), w\(r) = w(\r) (30) 

assumes the same boundary values at r = and r = oo, and the action S\ = S[m\, a\, w\] must have a 
critical point at A = 1, i.e. [dS/d\]\=i = 0. Therefore we find the following virial relation satisfied by the 
finite energy solutions of the field equations (note that following 33 j , it is possible to write a similar relation 
for black hole configurations, also) 

F : f * ' asgnsr lii)! w "^ d - *» - 1){2 " N ^ + <" - 2p - < 31) 



p=i J0 



12 2(d-3)m 2r 2 



+2p 72 -(^z^ + 1 t))=o. 
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For p = 1, i.e. a i ?2 -theory, the above relation reads 



dr 



o r^((d - 5)(W 2 + ^1!) + -' 2 (^(^ - 3) + £)) = 0, (32) 



which clearly shows that no nontrivial gravitating solution with finite mass exists for d > 4, since all terms 
in the integrant are strictly positive quantities. 

Therefore it becomes obvious that new terms in the YM-hierarchy should be introduced as the spacetime 
dimension increases. For a given d, the relation P > [(d + l)/4] should be satisfied. As it happens, to go 
to 5 < d < 9 it is necessary to include at least the second member of the YM hierarchy to provide the 
requisite scaling (similarly for 9 < d < 13 it is necessary to include the thrid member of the YM hierarchy). 
In practice we add only the lowest order such term necessary. 

We mention here also the Hawking temperature expression of the black hole solutions. For the line 
element Ipjl. if we treat t as complex, then its imaginary part is a coordinate for a nonsingular Euclidean 
submanifold iff it is periodic with period 

0= Nl( 4 w v (33) 
N'(r h )a{r h ) 

Then continuous Euclidean Green functions must have this period, so by standard arguments the Hawking 
temperature is (with ks = h = 1) 



T H = 



- (i _ S _^ + ( ,_ 1) £) s £_( i _ 1+(i _ 1) g). ,34) 



Thus the Hawking temperature of such systems appears to be suppressed relative to that of a vacuum black 
hole of equal horizon area. 

In the presence of higher oder terms in the YM action, dimcnsionlcss quantities arc obtained by rescaling 

r -> fa/n) 1 /* A — > (T1/T2) 1 / 2 , m(r)^ m (r) K (r 1 /r 2 )( <i - 3 )/ 4 1 (35) 

This reveals the existence of one fundamental parameter which gives the strength of the gravitational inter- 
action 

3/2 

a 2 = ^W, (36) 



1/2 ' 

KTr, 



and P — 2 independent coupling constants 



fe-i 

,2 T k I T\ ' 



Pt-2 = -[- . (37) 

with k = 3, P — 2 (i.e. no such constants appear in a p = 2 system). 

For the F(2) + F(A) systems in d = 6, 7, 8 considered in there exist gravity decoupling solutions at 
a — 0, from which the gravitating solutions branch out to a maximum value of a ma x and then decrease. 
The second limit of a — > also exists. For the same system in d = 5 studied in ^3] on the other hand, 
there exists no gravity decoupling limit but nonetheless the solution branches out from a = by employing 
a scaling procedure |19| . and again a increases to a value a max and decreases. But in this case it does not 
reach the a — > limit. Rather it stops at a new critical point a c around which it oscillates |14| . which is a 
new type of critical point identified in |19) and named a 'conical' fixed point. 

3.4 Numerical solutions 

In the present work, we restrict our attention to the simplest nontrivial cases with only the two terms p = 1 
and 2 in the YM hierarchy. However, we have obtained some numerical results also for a P — 3 hierarchy 



10 



with /?i = 1, which will be briefly mentioned. For d = 5, the solutions we found have some special features 
which will be discussed separately. 

Both regular and black hole solutions of the EYM-hierarchy equations appear to exist for any value of 
A. Given (a, d, A), AAdS solutions may exist for discrete set of shooting parameters b and Wh respectively. 
We follow the usual approach and, by using a standard ordinary differential equation solver, we evaluate the 
initial conditions at r = 10~ 5 (or + 1CP 5 ) for global tolerance 10~ 14 , adjusting for shooting parameters 
and integrating towards r — ► oo, and looking for AAdS solutions with a finite mass. 

Similar to the p = 1, d — 4 asymptotically flat case, it can be proven that for all AAdS (or A = 0) 
solutions, w(r) is confined within the strip < 1. This can proven as follows: we suppose the existence 

of solutions with w(r) > 1 for some interval of r. Therefore w must develop a maximum for some ro, 
w '( r o) — and io(rn) > 1 with w"(ro) < 0. However, the equations l|24[l imply that in the region w > 1 
the only extremum can be a minimum. Therefore the condition |xf(r) | < 1 is always fullfilcd. As a general 
feature, all solutions discussed in the rest of this section present only one node in the gauge function w(r). 
Similar to the A = case, we could not find multinode solutions 5 . 

The absence of multinode solutions in this F(2) + F(4) model with A = 0, in the relevant dimensions 
5 < d < 8 is analytically explained in ■ We expect that the relevant fixed point analysis yields qualitatively 
similar results also for |A| > 0. This is borne out by our numerical results. 

For any regular solution, the metric functions m(r) and cr(r) always increase monotonically with growing 
r from m(0) = and cr(0) = <To at the origin to m(oo) = M and ct(oo) = 1, respectively. The gauge 
function always interpolates between w(0) = 1 and w(oo) = — 1 without any local extrema. For black hole 
configurations, the behaviour of the functions m(r),a(r) and w(r) is similar to that for regular solutions. 
The gauge potential w(r) starts from some finite value < wir^) < 1 at the horizon and monotonically 
approaches —1 at infinity. The metric functions m(r) and cr(r) increase also monotonically with r. In the 
asymptotic region, the geometry corresponds to a Schwarzschild-AdS solution. However, although these 
solutions are static and have vanishing YM charges (w 2 (oo) = 1), they are different from the Schwarzschild- 
AdS black hole, and therefore are not fully characterised by the mass-parameter M. 

3.4.1 Regular solutions d = 5 

As a 2 — > 0, the YM equations present a nontrivial, finite energy solution in a fixed AdS background. This 
nongravitating configuration approaches in the A — limit the YM instanton in four dimensional flat space 

M- 

When a 2 increases, this solution gets deformed by gravity and the mass M decreases. At the same time, 
both the value cr(0) and the minimal value N m of the function N(r) decrease, as indicated in Figure 3. 
This branch of solutions exists up to a maximal value a max of the parameter a, which is smaller than the 
corresponding value in the asymptotically flat case For example, we find numerically a„ ax ~ 0.3445 
for A = —0.2 while the corresponding value for A = —0.01 is o^ ax ~ 0.5322. (without a cosmological term, 
this branch extends up to a^ ax « 0.5648.) 

Similar to the EYM theory with A = ^1], we found always another branch of solutions on the interval 
a 2 e [acr(i)' a max] wrtn a cr(i) depending again on the value of A (e.g. a 2 r ^ ss 0.2876 for A = —0.2). 
On this second branch of solutions, both <r(0) and N m continue to decrease but stay finite. However, a 
third branch of solutions exists for a 2 S a cr(2)] > on which the two quantities decrease further. A 

fourth branch of solutions has also been found, with a corresponding a 2 r ,^ close to ot 2 cr , 2 y Further branches 
of solutions, exhibiting more oscillations very likely exist but their study is a difficult numerical problem. 
Along this succession of branches, the main observation is that the value <r(0) decreases much faster than 
that of N m as illustrated in Figure 3. Also, the mass parameters do not increase significantly along these 
secondary branches. However, the shooting parameter b increases to very large values. The pattern strongly 
suggests that after a finite (or more likely infinite) number of oscilations of <r(0), the solution terminates into 
a singular solution with cr(0) = and a finite value of N(0). 

This is the behaviour observed in Jl] for the EYM theory with A = 0. The inclusion of a negative 
cosmological constant does not seem to qualitatively change the properties of the system, but leads to 

5 Multinode solutions exist if the lowest order YM term is F(2p) with p > 2, in the appropriate dimensions 1191 d. 
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different values of the critical parameters. As in the A = case, the dominant term at the gravity decoupling 
limit a — > is the F(2) term, the energy being given by the action of the usual instanton |39| . while as 
a — > a cr the dominant term is F(A). The mechanism for this effect is explained in |19) (for A = 0) and is 
supported by our numerical results here. The typical d = 5 globally regular solutions look very similar to 
the d = 6,8 profiles presented in Figures 4, 5. 

3.4.2 Regular solutions d > 5 

The solutions in this case resemble again the A = situation. In the presence of suitable higher order term in 
the hierarchy, the YM equations admit finite energy solutions in a fixed AdS^ A p — 1, AdS4 exact solution 
was found in 001 1 its higher order generalisation (specifically for p — 2, AdSs) for d > 4 being discussed in 
|41|. According to the standard arguments, this AdS soliton can be generalised in the presence of gravity, 
provided that the dimensionless coupling constant a is small enough. Therefore the gravitating solutions 
exist up to a maximal value a max of the gravitational coupling constant. This value a max for a given d 
depends on the value of A. For example in d = 6, a max (A = 0) = 0.12675; a max (A = —1) = 0.070422; in 8 
dimensions we find a max (A = 0) = 0.002193 while a max (A = —5) = 6.69 x 10~ 4 . 

When a increases, the mass of the gravitating solutions decreases while the function N(r) develops a 
local minimum N m which becomes deeper while gravity becomes stronger and the value cr(0) decreases from 
one. At the same time, the value of the shooting parameter b increases with a. Our numerical analysis for 
d < 10 indicates that a second branch of regular solutions always exists, starting at a max . Along this second 
branch the values <x(0) and N m decrease monotonically with a, while b and M still increase. The mass of a 
second branch solution is always larger than the corresponding mass (for the same value of a) on the first 
branch. For d > 6, the numerical analysis suggests that this second branch persists up to a 2 ~ and that 
in this limit cr(0) approaches a very small value. As far as our numerical analysis indicates, the value N m 
tends to a finite value in this limit so that there occurs no horizon. Therefore two regular solutions seems 
always to exist for any a < ad- 

The case d = 6 is special, since the numerical procedure fails to give reliable results for second branch 
solutions, starting with some ad, whose value is A-dependent (for example we found a^A = —1) = 0.0435 
while a<j(A = 0) = 0.0573). The quantity cr(0) reaches a very small value as a — > ad- The minimal value 
of N(r) remains finite so that no horizon is approached. We expect that a different parametrisation of the 
metric and variables would allow us to continue this second branch to a — > in this case, too. 

The behaviours just described qualitatively duplicate those of the A = case ^3], an d are analysed in 
|19j . Likewise, the solutions are dominated by the F(2) terms in the gravity decoupling limit a — > 0, while 
at the other end (on the second mass-branch) it is the F(4) terms that dominate. 

Typical d = 6, 8 solutions are presented in Figures 4, 5 respectively. In Figures 5, 6 and 7 we plot some 
relevant quantities for d = 6, 7, 8 and several values of A. One can see that the qualitative behaviour of the 
functions to, er, w does not change by changing the value of A. 

The results we found by including the p = 3 term in the YM-hierarchy for d = 9, 10 follows the same 
pattern. Although the picture gets more complicated by the existence of one more coupling parameter, two 
branches of solution are always found to exist. We noticed also the existence of a maximal value of a which 
is (A,/3i) dependent. 

3.4.3 Black hole solutions d = 5 

According to the standard arguments, one can expect to find black hole generalisations for any regular 
configurations, at least for small values of the horizon radius rv For completeness we discuss here the basic 
features of the AAdS black hole solutions. 

Again, the case d = 5 is special. The properties of these AAdS solutions are rather similar to the five 
dimensional asymptotically flat black hole configurations discussed in ^3]. First, black hole solutions seem 
to exist for all values of a for which regular solutions could be constructed. Also, solutions exist only for a 
limited region of the (77, , a) space. 

The typical behaviour of solutions as a function of is presented in Figure 9, for a small value of a 
compared to the maximal value a max of the regular solutions, in which case we notice the existence of only 
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one Th = regular configuration. Starting from this regular solution and increasing the event horizon radius, 
we find a first branch of solutions which extends to a maximal value rh( max y As seen in Figures 9, 10, 
the value of rhi ma x) depends on A, a. The Hawking temperature decreases on this branch, while the mass 
parameter increases; however, the variation of mass and cr{rh) is relatively small. Extending backwards in r^, 
we find a second branch of solutions for < r h ^ max y This second branch stops at some critical value r h ^ cr y 
where the numerical iteration fails to converge. The value of cr(r^) on this branch decreases drastically, as 
shown in Figure 9. Also, the Hawking temperature after initially increasing, strongly decreases for values 
near r/j( cr ), approaching a very small value, while the increase of the total mass is still very small. Similar 
to the A = case [T3|, higher branches of solutions on which the value er(r/,) continues to decrease further 
to zero are likely to exist. However, the extension of these branches in r/, will be very small, which makes 
their study difficult. An approach to this problem with different parametrisation appears to be necessary. 

However we find, that the global picture is changed by considering large enough values of a. In this case, 
more than one regular configuration exists for a given value of a. This situation is illustrated in Figure 10, 
for solutions with a 2 = 0.5. Two regular solutions exist for this particular value of a and we find two black 
hole branches connecting these = configurations. Again, the mass of the second branch solutions is 
always larger than the corresponding mass on the first branch. 

Preliminary numerical results indicate an even more complicated picture for solutions with a near a max . 
In this case the configurations combine features of both types of solutions discussed above. Several branches 
of black hole solutions are found for the same a. These branches start from regular configurations and are 
possibly disconnected. 

3.4.4 Black hole solutions d > 5 

Although predicted in no discussion of the A = 0, d > 5 black holes is presented in literature. 

Again, black hole counterparts appear to exist for any regular solution. However, solutions with the right 
asymptotics are found for a limited region of the (r^, a) space only. We plot in Figures 11-13 some results we 
found for d — 6, 7 and d = 8. Starting for a given ao < a max from a = first branch regular solution, we 
found the existence of a branch of black hole solutions extending up to a maximal value of the event horizon 
radius r% = r™ ax . When increases, the mass and the Hawking temperature increase while the value u(rh) 
decreases from its value at the regular solution. A second branch of black hole solutions seems to appear 
always at r™ ax , extending backwards to a zero event horizon radius. This limiting solution corresponds to 
the second branch of the regular solution at this value of a = ao. Like in those d = 5 cases where there 
are two regular solutions for a given a, say ao, here too these two regular solutions are the ru — > limits 
of the corresponding black hole "loop", connecting the two regular solutions. Along this second branch 
the values cr(r/ l ) of the metric function a on the event horizon decrease monotonically with rn, while the 
Hawking temperature strongly increases. The mass of the solution of the second branch is larger than the 
corresponding one on the first branch, as illustrated by Figures 10-12. It is interesting to note here from 
the point of view of numerics, that a black hole loop corresponding to the two 77, — ► limits of two regular 
solutions, say at ao, can be constructed numerically even when the numerical process for the higher mass 
branch of the regular solutions for this ao runs into difficulties. 

The introduction of a negative cosmological constant does not appear to change this picture qualitatively. 
However, we notice a smaller value of r™ ax with increasing |A| and a larger value of M for the same a. 

In Figure 14, we present the profiles of the metric functions N and a and gauge function w(r) for the 
same values of a 2 , A on the first and second branch for some d = 6 solutions. The dependence of these 
functions on the value of A is illustrated in Figure 15 for several d = 7 black hole solutions. 

4 A computation of mass and action 
4.1 The counterterm method 

It is well known that the generalisation of Komar's formula for AAdS spacetimes is not straightforward and 
requires the further subtraction of a background configuration in order to render a finite result for the mass. 
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This problem was addressed for the first time in the eighties, with different approaches (see for instance 
Ref. Another formalism to define conserved charges in AAdS spacetimes was proposed in 01] and 

uses conformal techniques to construct conserved quantities yielding the results obtained by Hamiltonian 
methods. Other more recent approaches to the same problem are presented in Ref. |45| . 

As expected, these different methods yield the same total mass for the spherically symmetric AAdS 
configurations considered in Section 3 

AW = {d - 2) " d - 2 M, (38) 

S7TG 

where Q d _ 2 = 2TT < - d - 1 ^ 2 /T({d - l)/2) is the area of a unit (d — 2)-dimcnsional sphere, and M is defined in 
the second member of (|29[) . 

A procedure leading (for odd dimensions) to a different result has been proposed by Balasubramanian 
and Kraus |2J- This technique was inspired by AdS/CFT correspondence and consists in adding suitable 
counterterms I c t to the action of the theory in order to ensure the finiteness of the boundary stress tensor 
derived by the quasilocal energy definition |46j . These counterterms are built up with curvature invariants 
of a boundary dAi (which is sent to infinity after the integration) and thus obviously they do not alter the 
bulk equations of motion. Unlike background subtraction, the counterterm approach does no require the 
identification of a reference spacetime. Given the potential relevance of the EYM solutions in an AdS/CFT 
context, we present here a computation of the boundary stress tensor and of the total mass, by using the 
counterterm prescription. 

The following counterterms are sufficient to cancel divergences in a pure gravity theory for d < 9, with 
several exceptions (see e. g. [17] ) 



8nG 



d d - x xy 
8M 



d-2 £6(d-4) £ 3 ®(d~6) / AB d-1 



-K 



( R " K ' ,B -3& 2 ) < 39 > 



2(<2-3) 2(d- 5)(d- 3) 2 V ™ 4(d-2 
£ 5 Q (d - 8) ( U + 2 _ nn AB UAB _ d(d + 2) ^ 



(d-2) 3 (d-4)(d-6) \4(d-l) 16(d-l) 

-2k ab k cd k ABC d - —^—v A nv A n + v c n AB v c n AB 

i(d- 1) 

Here TZabcd, T^-ab are the Riemann and Ricci tensors, 1Z is the Ricci scalar for the boundary metric h and 
Q(x) is the step function, which is equal to 1 for x > and zero otherwise; A, B,.., indicate the intrinsic 
coordinates of the boundary. 

Using these counterterms one can construct a divergenceless boundary stress tensor, which is given by 
the variation of the total action at the boundary with respect to h ab- Its explicit expression, restricting for 
simplicity to d < 7, is 

Tab = ^(K A b - Kh AB - ^^h AB + -—E AB ), (40) 
87rG I d — 3 

where Kab = — An B + V B n A ) is the extrinsic curvature defined in terms of the normal i%a to the 
boundary, K is its trace, and E AB is the Einstein tensor of the intrinsic metric Ha B - (The corresponding 
form of lHUI) for d = 7, 8 is given e.g. in (15]1. 

The result we find in this way for T AB is given by 



rpB _ 

L A — 



-i -p/ 2p— 1 \ 



where ua = S A and p is an integer. We can use this approach to assign a mass-energy to an AAdS geometry 
by writing the boundary metric in an ADM form 

h AB dx A dx B = -N^dt 2 + (7 ab (dx a + N^dt){dx h + N b a dt) (42) 
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and the definition of the energy in this context is 



E = [ d^xy^Nze. (43) 

Here e = u^u^T^ is the proper energy density while u M is a timelike unit normal to S. 

If there are matter fields on M., additional counterterms may be needed to regulate the action. This is 
the case of F(2) theory, discussed in Appendix B. The couinterterm action depends in this case not only on 
the boundary metric but also on the boundary value of the gauge field. 

However, we find that for P > 1 EYM solutions with A < in d = 5, 6, 7, 8 dimensions, the prescription 
(I39[l removes all divergences. The use of higher order terms in the YM curvature, namely F(2p) forms with 
p > 1 introduced in Section 3, results in this regularising of the masses. The crucial point here is that these 
solutions approaches asymptotically a Schwarzschild-AdS background, and the YM asymptotic parameter 
Wi appears only in the next to leading order of the Tab expression. 

The mass-energy of solutions computed in this way is 

E= {d ~SnG~ 2 M + E ° (44) 



where, for 3 < d < 9 



The additional term Eo appearing in E for d — 5, (7) is the mass of pure global AdS^^-j) and is usually 
interpreted as the energy dual to the Casimir energy of the CFT defined on a four (six) dimensional Einstein 
universe |21| . 

The metric restricted to the boundary \iab diverges due to an infinite conformal factor r 2 j£ 2 . The 
background metric upon which the dual field theory resides is 

I 2 

JAB = lim -irflAB- (46) 

For the asymptotically AdS^ solutions considered here, the (d — 1) dimensional boundary is the Einstein 
universe, with the line element 

lAB dx A dx B = -dt 2 + fdn 2 d _ 2 . (47) 

In light of the AdS/CFT correspondence, Balasubramanian and Kraus have interpreted Eq. I|40(l as 
< t ab >= -j= , where < t ab > is the expectation value of the CFT stress tensor. Then, the 
divergences which appear are simply the standard ultraviolet divergences of a quantum field theory and we 
can cancel them by adding local counterterms to the action. Corresponding to the boundary metric l|47|) . 
the stress-energy tensor tab for the dual theory can be calculated using the following relation 0^] 

.AB ^ ^_ i- I TuAB 



< t BC >= lim y —hh Tbc ■ (48) 



4.2 Action and entropy 

The above approach can be used to compute the Euclidean action and to prove in a rigorous way that the 
entropy of the EYM black hole solutions is one quarter of the event horizon area, as expected. Here we start 
by constructing the path integral ^Hj 

Z = [ D[g]D[V]e- u teW (49) 
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by integrating over all metrics and matter fields between some given initial and final hypersurfaces, \& 
corresponding here to the SU(2) potentials. By analytically continuing the time coordinate t — ► ir, the path 
integral formally converges, and in the leading order one obtains 

Z ~ e- 1 " 1 (50) 

where I c t is the classical action evaluated on the equations of motion of the gravity/matter system. The 
physical interpretation of this formalism is that the class of regular stationary metrics forms an ensemble of 
thermodynamic systems at equilibrium temperature T (see e.g. Z has the interpretation of partition 

function and we can define the free energy of the system F = — /3 _1 log Z. Therefore 

log Z = -I3F = S -0E, (51) 

or 

S = f3E-I cl , (52) 

straightforwardly follows. 

To compute I c i, we make use of the Einstein equations, replacing the R — 2 A volume term with 2R\ — 
16t:GT}. For our purely magnetic ansatz, the term T\ exactly cancels the matter field Lagrangean in the 
bulk action. The divergent contribution given by the surface integral term at infinity in R\ is also canceled 
by /surface + let and for 3 < d < 9 we arrive at the simple finite expression 

Replacing I c i now in l|52|) (where E is the mass-energy computed in Section 4.1), we find 

S=\n d _ 2 rt 2 , (54) 

which is one quarter of the event horizon area, as expected. 

^From the AdS/CFT correspondence, we expect the nonabelian hairy black holes to be described by some 
thermal states in a dual theory formulated in a Einstein universe background. The spherically symmetric 
solitons will correspond to zero-temperature states in the same theory. The existence of these hairy config- 
urations suggest that there should be some observables in the dual CFT that encode the hair information. 



5 Conclusions 

Motivated by recent results in EYM theory in four dimensional AAdS spacetime, we studied higher dimen- 
sional spherically symmetric solutions with non Abelian fields in the presence of a negative cosmological 
constant. Since the mass-energy of the AAdS configurations plays a central role in its application to the 
AdS/CFT correspondence, we emphasised this aspect of the solutions we found. The mass-energy of the 
usual EYM solutions (both asymptotically flat and AAdS), defined according to the standard prescription, 
always diverges in spacetime dimensions d > 4. One of the tasks performed in this work was a demonstration 
of this fact in the AAdS case. The main properties of these higher dimensional F(2) solutions resemble the 
d = 4 case, a continuum of solutions with arbitrary asymptotic values of the gauge function w(r) being 
found. 

Then we focused on two possible approaches of dealing with the divergences of the mass and action. One 
of these involved the regularisation of the mass-energy using a counterterm mechanisms. In this approach, it 
turns out that the counterterm action depends not only on the boundary metric, but also on the boundary 
values of the gauge fields. However, the masses of the solutions defined in this way may take negative values, 
leading us to the second approach. 
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The other, which formed the main thrust of the work, was to augment the action density of the system 
with higher order curvature terms, consisting of 2p-form curvatures F(2p). These terms were added to the 
usual YM system constructed from F(2). It resulted in EYM solutions supporting finite mass-energy in all 
spacetime dimensions 5 < d < 2p. 

Concerning the construction of regular finite energy classical AAdS solutions in higher dimensions, we 
restricted ourselves to systems consisting exclusively of gravitational and non Abelian gauge fields. The 
salient features of the resulting solutions are captured in this framework, the addition of other (string theory 
inspired) matter terms being deferred to later work. The asymptotically flat versions of the higher-p EYM 
systems having been studied in ^J^J^j, our task here involved the introduction of a negative cosmological 
constant. The most succint way of listing our conclusions is: 

• The qualitative properties of the regular AAdS solutions in spacetime dimensions d = 5, 6, 7, 8 are the 
same as the corresponding asymptotically flat ones, namely 

— A one parameter family of solutions parametrised by the (dimensionless) gravitational coupling 
constant a start at a = (the gravity decoupling limit) and exist up to a maximum a ma x) after 
which a decreases again and ends at critical value. 

— For d = 6, 7, 8 the value of a in the second, not that of gravity decoupling, endpoint becomes very 
small and stops. In the d = 8 case it actually vanishes. For d — 5 the value of a in the second 
endpoint reaches a finite critical value, where it does not stop, but oscillates around this critical 
value. 

— As long as the physically important F{2) term is present in the YM sector, there exist no multinodc 
solutions. 

• Black hole counterparts appear to exist for any regular solution. The qualitative properties of the 
AAdS black hole solutions are similar to the asymptotically flat case: 

— Different from the four dimensional theory, the event horizon radius presents a maximal value. 
This maximal value is a function of the gravitational coupling constant a. 

— For d = 6, 7, 8 the black hole solutions form a loop connecting the two regular solutions with the 
same value of a. The solutions of the five dimensional theory are somehow special, presenting a 
complicated brach structure which depends on a. 

Axially symmetric generalizations of these solutions are likely to exist. We expect them to share the 
basic properties of the spherically symmetric configurations. 
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Appendix A: A nonexistence proof for p = 1, d > 5 finite mass 
solutions with A < 

Following the notation used in we introduce a new variable 

z = 21ogr, (55) 
and rewrite the basic equations (JBJl, © in the form 

dm Ti (d-5) . / dw. n (d — 3) . 2 ,n2 



dz 4 e^^(_)^^K-irj, (56) 

(d — 5)N (.d-3) _m e z T\e~ z 2 2 \ (d — 3) . 2 



where the metric function a has been eliminated by using i|10|) . The function N is given by JfjJ), namely as 
TV = 1 - |e-^ 3 ) z / 2 m(z) + e z /e. 

To devise a proof for the nonexistence of finite mass solutions of the above system, it is convenient to 
introduce the function 

which, from l|56|) and (|57|l satisfies the equation 

dE ( dw\ 2 ( Vi _ Z1 ^ e z m,, „. (d-3) \ 

The relations in Appendix B of Ref. are recovered for d = 5, where a nonexistence proof is presented 
for particle-like solutions with finite mass (the extension to black hole case is considered in [HI])- Therefore 
in what follows we will take d > 5 only. 

The approximate form of the function E at the origin (or event horizon) and infinity is found by taking 
P = 1 in relations l|26l) . (|27|l and (|29|l given in Section 3 and replacing in (|58|) . It is obvious that E{rh) < 0, 
since N(rh) = 0; at the same time, the corresponding expression as r — > (i.e. z — * — oo) is 

E~ ^(5-d)b 2 e 2z + (60) 

and we find E — > — in this limit. 

Also, the relations Q26)). I|27|) give m(r = 0) > 0, m(rh) > 0, which, together with the equation (|56|l 
implies that the mass function m(r) is positive definite. 

Besides, by replacing in (|58|l the asymptotic expressions (|29|) as r — > oo it follows that 

E ^ ^^ e _^ z + l_ wl[d m 5)e _ {d _ 3)z + (6i) 

i.e. E — > +0 as z — > oo, for finite mass solutions. 

Therefore, if the solution is regular everywhere, E must vanish at some finite point zq, and dE / dz > 
there, with E > for z > z$ (when there are several positive roots of E, we take the largest one). However, 
another point should exist z\ > zq such that dE/dz = i.e. the function E should present a positive 
maximum for some value of z. Now we integrate the equation (|59|l between zq and z\ and find 

' / dw\ 2 ( t\ e z m . (d-3) . 



EM = ~J ro [^) ^e- z E + - + -(d^)e-^ z )dz<0, (62) 

which contradicts E(z\) > 0. Therefore E(z) should vanish identically and one finds no d > 5 finite mass, 
spherically symmetric EYM configurations in a F(2) theory. Note that this argument does not exclude the 
existence of configuration with a diverging mass functions as r — > oo. 
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Appendix B: A matter counterterm proposal 



In this Section we comment on the issue of mass definition of AAdS solutions in a F(2) 2 (i.e. p = 1) 
theory, if we do not exercise the option of employing higher order YM survature terms. As proven in Section 
2, although the spacetime is still AAdS, the mass function m(r) of these solutions generically diverges as 
r d-5 ^ or ag i g r f or d = 5). AAdS solutions with a diverging ADM mass have been considered recently by 
some authors, mainly for a scalar field in the bulk (see e.g. I n this case it might be possible to relax 

the standard asymptotic conditions without loosing the original symmetries, but modifying the charges in 
order to take into account the presence of matter fields. 

Similar to the case of scalar field, for d > 5 it is still possible to obtain a finite mass of EYM solutions 
in a F 2 theory by allowing I ct to depend not only on the boundary metric \iab-, but also on the gauge field 
strength tensor. This means that the quasilocal stress-energy tensor (|40J) also acquires a contribution coming 
from the matter fields. 

We find that by adding a counterterm of the form 

= -j^— ( d^xV^h tr F AB F AB (63) 
yd - 5) J aM 

to the expression i|39|) . the divergence disappears. This yields a supplementary contribution to l|40|) 

Tab = - _ * , 1 , hAB tr F CD F CD . (64) 
07rG d — 5 

The mass of the d > 5 solutions computed in this way is 

E = {d - 2) ° d - 2 M + E (65) 

07T(jr 

where Mq is the constant appearing in the asymptotic expansion (|16fl . It can also be proven that this 
prescription leads to a finite action and the entropy-area relation is satisfied. However, as seen in Figure 2, 
the parameter Mq of the p = 1 solutions takes negative values, pointing to some pathological properties. 

It would be nice to have a rigorous derivation of the matter counterterm expression, possibly along the 
lines of Ref. Also, there remains the issue of the d = 5 solutions, whose logarithmic divergeces require 
a different approach. 
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Figure 2. The value N m of the minimum of the metric function N(r), the parameter M$ appearing in 
the asymptotic expansion of m(r), the asymptotic value wo of the gauge function w(r) as well as the value 
cr(0) of the metric function a at the origin, are shown as functions of b for d — 6 solutions of F 2 theory with 
A = -1. 
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Figure 3. The value N m of the minimum of the metric function N(r), the mass parameter M as well 
as the value of the metric function a at the origin, er(0), are shown for d — 5 solutions as functions of a 2 
and several values of A. 
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Figure 4. The profiles of the functions m(r), cr(r) and w(r) are plotted as functions of radius for typi- 
cal d = 6 regular solutions in a EYM theory with p = 1,2 terms and a 2 = 0.1, A = —0.01. Here and in 
Figure 14 the continuos line corresponds to an upper branch solution, the dotted line denoting lower branch 
profiles. 
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Figure 5. Typical globally regular solutions of d = 8 EYM theory with a 2 — 0.0005 and p = 1, 2 terms are 
plotted as functions of radius for several values of A. 
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Figure 7a. 
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Figure 8a. 
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Figure 9a. 
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Figure 9b. 

Figure 9 In Figure 9a we plot the value of the gauge field function at the horizon w{rh) and a(rh), the 
value of the metric function a at the horizon (the magnified profiles of cr(rh) are displayed in the window, 
to help distinguish these from the profiles of w(r^)). In Figure 9b, the mass parameter M and the Hawking 
temperature Th are presented. All profiles as functions of the event horizon radius for the p = 1, 2 black 
hole solutions in five dimensions with a 2 = 0.02 and several values of A. 
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Figure 10a. 
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Figure 10b. 

Figure 10 The value of the gauge field function at the horizon w{rh) and d(r^), the value of the metric 
function a at the horizon, (Figure 10a), as well as the mass parameter M and the Hawking temperature Th 
(Figure 10b) are shown as functions of the event horizon radius for the p = 1, 2 black hole solutions in 
five dimensions with a 2 = 0.5 and several values of A. 
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Figure 11. The value of the gauge field function at the horizon w(rh) and a(rh), the value of the met- 
ric function o at the horizon, as well as the mass parameter M and the Hawking temperature Th are shown 
as functions of the event horizon radius ru for the p = 1, 2 black hole solutions in six dimensions with 
a 2 = 0.066 and several values of A. 
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Figure 12. The value of the gauge field function at the horizon w(rh) and a(rh), the value of the met- 
ric function o at the horizon, as well as the mass parameter M and the Hawking temperature Th are shown 
as functions of the event horizon radius r/j for the p = 1, 2 black hole solutions in seven dimensions with 
a 2 = 0.011 and several values of A. 
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Figure 13. The value of the gauge field function at the horizon w(rh) and a(rh), the value of the met- 
ric function a at the horizon, as well as the mass parameter M are shown as functions of the event horizon 
radius for the p = 1, 2 black hole solutions in eight dimensions with a 2 = 0.002 and several values of A. 
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